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Striped phases, in which spin, charge, and pairing correlations vary inhomogeneously in the Cu02 
planes, are a known experimental feature of cuprate superconductors, and are also found in a 
variety of numerical treatments of the two dimensional Hubbard Hamiltonian. In this paper we use 
determinant Quantum Monte Carlo to show that if a stripe density pattern is imposed on the model, 
the d-wave pairing vertex is significantly enhanced. We attribute this enhancement to an increase in 
antiferromagnetic order which is caused by the appearance of more nearly half-filled regions when 
the doped holes are confined to the stripes. We also observe a 7r-phase shift in the magnetic order. 

PACS numbers: 71.10.Fd, 02.70.Uu 



I. INTRODUCTION 

Cooper pairs in conventional BCS superconductors are 
typically envisioned to have a large spatial extent char- 
acterized by the coherence length ^, which is many hun- 
dreds of lattice spacings in elemental, metallic supercon- 
ductors. At the other extreme is the BEC regime, where 
much smaller Cooper pairs form into bosonic particles 
which subsequently can Bose-Einstein condense into a 
superfluid phase. The crossover between the BCS and 
BEC limits has been much explored^'. 

If real space pairing on small length scales is fa- 
vorable energetically, a natural question to ask is why 
more than two fermions do not bind together. Indeed, 
the competition of such "phase separation" with su- 
perconductivity was a central theme of investigation in 
early studies of Hamiltonians like the two dimensional 
fermion HubbarcP'^ and t-fi^^^ models in the context 
of cuprate superconductors, which have short coherence 
length. Upon doping away from one particle per site, 
a spatial division into half-filled AF regions and areas 
where the hole concentration is high occurs. The ten- 
dency to phase separation was found to be especially 
great in the t-J model, and somewhat less so in the Hub- 
bard model. Roughly speaking, for small hole doping 
in the strong coupling limit such phase separation can 
be envisioned to arise because it minimizes the energy 
by preserving the largest number of antiferromagnetic 
bonds. On the other hand, when t » J and the amount 
of holes is large, phase separation was shown to occur as 
to minimize the doped hole kinetic energjA^i 

A compromise between complete phase separation into 
distinct two dimensional regions, and spatial homogene- 
ity which would be favored by longer range Coulomb 
interactions is the possibility that the doped par- 
ticles form quasi-one dimensional patterns in which hole- 
rich and hole-poor regions alternate. Magnetic domain 
lines were observed in inhomogeneous Hartree-Fock stud- 
ies very early in the history of cuprate superconductivity. 
'Charge domain lines' were first reported in a multi-band 



model which included both the copper d and oxygen p 
orbitals^^, and subsequently in the single band Hubbard 
Hamiltonian-2^. 

Such strip ed p atterns have also been observed 
experiment ally23l26|^ jjj Lai.6-a;Ndo.4Sra:Cu04 there is a 
suppression of superconductivity at x = 1/8 associated 
with a tilt pattern of oxygen tetrahedra and charge/spin 
domain walls. Away from a; = 1/8 domain wall ordering 
is weaker, with coexisting superconductivity. 

The relation between pairing and stripes is still con- 
troversial. The formation of charge stripe-order is re- 
lated to an increase in the resistance perpendicular to 
the CuO planes, frustrating the formation of bulk super- 
conductivity, as denoted by a sharp decrease of the crit- 
ical temperatur^^, Tc, resulting in values as low as 4 K. 
As stripes are perpendicularly stacked in each successive 
CuO plane, if the Josephson coupling between stripes in 
the same plane is negative thejxmpling between stripes 
in different planes is destroyecPsEU. 

In LaBaCuO, on the other hand, angle-resolved spec- 
troscopy and transport measurements sug gest a positive 
correlation between 2D pairing and stripegS^EIl. 

There is an order of magnitude dropP^'^^Ml [j^ 
in-plane resistivity, Pab, when the spin order occurs at 
Tso ^ 42 K. Furthermore, there are indications of true 
2D supe rcond uctivity for T < Tbkt ~ 16 K when pab 
goes to Cp332l. 

A recent exciting development is the direct evidence of 
stripes in YBa2Cu30y . Nuclear magnetic resonance mea- 
surements show that high magnetic fields induce charge 
order with the same period of four lattice spacings as 
in LaBa b ased cupratea^. Zero-field diffraction with res- 
onant softP^'^ and harcP^ x-ray scattering, observed in- 
commensurate charge order with a period around 3.2 lat- 
tice spacings. These experiments seem to indicate that 
stripes are an intrinsic phenomena in the cuprates, and 
that at least incipient charge ordering can be seen even 
in compounds with larger critical temperatures. 

This has led to a large set of theoretical studies of 
charge ordered patterns to refine and improve upon the 
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initial mean field treatmenlP^^. Given the complexity 
of the question of superconductivity even in the homoge- 
neous model, it is not surprising that the details of the 
interplay with stripe formation should be challenging. 

A recent calculation within the dynamical cluster ap- 
proximation (DCA)'^'' provided detailed information con- 
cerning pairing correlations amidst static stripes, and re- 
vealed a rich competition between an enhancement of 
the pairing interaction and a suppression of the non- 
interacting susceptibility. The latter effect occurs as a 
consequence of the formation of Mott regions away from 
the stripes, and a resulting suppression of quasiparticle 
weight. The two effects combine to give a non-monotonic 
evolution of the transition temperature with the stripe 
modulation strength. 

In this paper we undertake determinant Quantum 
Monte Carlo (DQMC) studies of stripe formation in the 
two dimensional Hubbard model. The DQMC method 
complements the DCA approach, working in real space 
rather than momentum space. It is possible to study 
somewhat larger clusters with DQMC, with, however, an 
off-setting greater restriction to the accessible tempera- 
tures. The key results of our work are an enhancement 
of antiferromagnetic order by charge modulation, and a 
TT-phase shift above a critical threshold of the stripe po- 
tential. Accompanying this larger magnetic order is a 
stronger signal of d-wave pairing in the associated super- 
conducting vertex, although we are not able to lower the 
temperature enough to cross below the transition tem- 
perature. 

II. HAMILTONIAN AND METHODOLOGY 

We study a two dimensional repulsive Hubbard Hamil- 
tonian in which stripes are introduced externally via a 
raised site energy Vq on a set of rows of period V, 

^ = E (^L^j. + 4-^-) + ^ E 
- M E(^^it + ^ii) + ^0 E ^^'^ + ^^i^) 

i iyEV 

Here i = 1 is the intersite fermion hopping between near 
neighbor sites i, j on a square lattice, U is an onsite repul- 
sion, is a global chemical potential, and Vq is an addi- 
tional on-site energy imposed on a set of rows i = {ix,iy) 
with inod{iy,V) = 0. When P = 4 this choice produces 
the spin and charge patterns postulated based on neutron 
scattering studies^^ and shown to arise in density ma- 
trix renormalization group (DMRG) studies on the t-J 
HamiltoniarPI'. 

This Hamiltonian does not, of course, address the is- 
sue of spontaneous stripe formation in a translationally 
invariant system, nor does it acknowledge the tendency of 
charge domain walls to fluctuate. Nevertheless it allows 
us to examine the nature of spin and pairing correlations 
in the presence of a set of pre-formed lines of reduced 
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FIG. 1: (Color online) 16x16 lattice with period "P = 4, sites 
with Vo active are depicted in blue (filled cricles), whereas the 
interstripe sites are red (empty circles). 



charge density, and is an appropriate approximate model 
in the limit where the energy scale of stripe formation is 
considerably greater than that of pairing. 

Most of our results will be for 16x16 lattices with V = 
4, so that each stripe (row with Vq term active, blue filled 
circles in Fig. [T]) is separated by three rows where the Vq 
term is not present (empty circles in Fig. [T]) . The entire 
16x16 site system accomodates four stripes for this V = 4: 
case. We will also analyze finite size effects and present a 
smaller amount of data for pairing correlations at other 
periodicities V. Most of our results will be for a total 
density p = 0.774 averaged over the entire lattice. This 
value was chosen because it allows for the exploration 
of a broad range of densities on the stripe and between 
stripes, and also because certain experimentally observed 
characteristics of the striped phase, such as the 'tt phase 
shift' of the spin correlations when traversing a stripe, 
are absent at some other fillings. 

In this "P = 4 case, we expect the charge order to 
be modulated along the y direction with the same pe- 
riod four, and the spin order has a period twice as large 
as a result of the '7r-phase shift' (see below), also ob- 
served in experiments. Therefore, X-ray and neutron 
diffraction experiments are expected to show peaks at 
Qco = ^(0;±2,S) and Qso = ^{0.5;0.5 ± 5^^^, 
respectively. The value of 6 is doping dependent in 
the La2-2:Ba2,Cu04 familj^^^and is assumed to be 1/8 
for doping higher than 0.126^^^ in agreement with the 
stripe sketch in Fig. [T] Some experiments, however, 
displayed incommensurate charge stripes with S smaller 
than expectecP^*^. 

Our methodology is determinant Quantum Monte 
Carlo (DQMC)'"^^ In this approach, the quartic in- 
teraction term is replaced by a coupling of the local z 
component of spin to an auxiliary field^^. The fermion 
degrees of freedom are integrated out analytically, leav- 
ing a Monte Carlo over the auxiliary field. In the process 
of eliminating the interaction term the inverse tempera- 
ture (3 is discretized. We choose the discretization mesh 
At = l/8t. The resulting Trotter errors are typically 
only a few percent, and their elimination is of conse- 
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quence only if subtle changes in short range correlation 
functions are of interest^® In DQMC, the Trotter er- 
rors are typically smaller than accessible statistical errors 
on long range spin, charge and pairing correlations at the 
lowest temperatures. 

The stripe potential Vo breaksparticle-hole symmetry 
so that there is a sign problenPSI for all fillings. Thus 
we cannot obtain ground state properties as can be ac- 
cessed, for example, in the half-filled homogeneous sys- 
tem or the attractive Hubbard Hamiltonian at any filling. 
A contour plot of the sign for p = 0.774 and [7 = 4 for 
a 16x16 lattice, is shown in Fig. [2] It can be seen that 
the sign problem is no worse than in the doped, homo- 
geneous case where one can, already discern significant 
trends concerning the superconducting correlations^^. 

We will show results for the spin and (d-wave) pair 
correlations, 
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The quantities Cspin(i) and Cdpair(i) do not have com- 
plete translation invariance owing to the presence of the 
stripes. They will depend on the row index iy, for ex- 
ample whether mod(zy,'P) = so the row has a reduced 
density, or, for mod(iy,'P) ^ on the distance from the 
reduced density stripe. 

In a phase with long range spin or pairing order, the ap- 
propriate correlation function would approach a constant 
value asymptotically as |i| — ?► oo. Indeed, precisely this is 
seen in the ground state of the half-filled Hubbard model 
which has long range anti-ferromagnetic ordei^^'^, and 
the attractive Hubbard model at a range of filling^^. 
It is still an open issue whether Crfpair(i) is nonzero at 
long distances when the homogeneous Hubbard model is 
doped away from half-filling, because the sign problem 
prevents attaining the ground state in the doped case. 

Given this limitation on the simulations, it is impor- 
tant to develop methods which extract the maximal use- 
ful information about the tendency to order at temper- 
atures above the putative superconducting phase tran- 
sition. To this end, one introduces a somewhat more 
sensitive measure of pairing by considering the pair-field 
susceptibility Pd and its associated vertex. To define Pd 
we first extend the definition of the equal time pair cor- 
relation function Cdpair(i) to allow the insertion and re- 
moval of the Cooper pair to be separated in imaginary 
time. Pd is the sum over all spatial sites i and integral 
over all imaginary time t separations of Cdpair(i, t). 



pair (^i ^) 



Pd = 



(A,j+i(r)At.(0)) 

y2 / Cdpair(i, t) dr 
i Jo 



(3) 



FIG. 2: (Color online) Contour plot of the sign as a function 
of Vo and /3 for a 16x16 lattice with f/ = 4 and p = 0.774. 
The contours are basically vertical, indicating that the sign 
problem is independent of Vo- However for Vb ~ 4 there is 
a modest improvement in the sign. This is the potential for 
which the interstripe density goes through half-fiUing. 



We also define the uncorrelated pair field susceptibil- 
ity Pd which instead computes the expectation values 
of pairs of operators prior to taking the product, with 
expressions like ( Ci_|_j^(r) C;^j^(T) Cjj.(0) Cj^(O) ) which 
appear in evaluating the Pd in Eq. [3| being replaced by 

Pd includes both the renormalization of the propaga- 
tion of the individual fermions as well as the interaction 
vertex between them, whereas Pd includes only the for- 
mer effect. In short, in DQMC, the averaging over the 
Hubbard-Stratonovich field replaces the interaction with 
the one body potential by the original electron-electron 
interactions, so that the order of averaging and multi- 
plying the operators can be used to control which many 
body effects are included. 

By evaluating both Pd and Pd we are able to extracl!^ 
the interaction vertex Td, 



Td = 



1 

Pd 



1 

7d 



(4) 



If TdPd < 0, the associated pairing interaction is attrac- 
tive. In fact, rewriting Eq. |4]as. 



Pd 



Pd 



1 + TdPd 



(5) 



suggests that TdPd — 1 signals a superconducting in- 
stability. This is the analog of the familiar Stoner crite- 
rion Uxo — 1 which arises from the random phase ap- 
proximation expression X — Xo/(l~C^Xo) for the inter- 
acting magnetic susceptibility x in terms of the noninter- 
acting xo- We will discuss this criterion in more detail in 
the coming sections. 
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FIG. 3: (Color online) The density of particles is shown 
on the striped rows where Vb acts and the unstriped rows. 
The total density of the system is fixed at p = 0.774 (circles 
for 16x16 and up triangles for 8x8 lattices), and p = 0.875 
(squares for 16x16 and down triangles for 8x8 lattices), There 
is only a very small variation of density on the unstriped rows 
with different distance from the stripes, so only the average 
is shown. Data for 8x8 and 16x16 lattices are essentially in- 
distinguishable. 
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FIG. 4: (Color online) The spin correlation function Cspin(i) 
down the center of one of the (three site wide) interstripe 
regions. Here i7 = 4, and the average density p = 0.774 
over the whole lattice, for Vb = 4 (a) and Vb = 10 (b). The 
interstripe region is fairly close to half-filling, and so, as the 
temperature T is lowered, Cspin(i) oscillates over fairly large 
distances. The arrows indicate the sites along which the spin- 
spin correlations are calculated. 



III. RESULTS 

When the total density, averaged over the entire lat- 
tice, is fixed, the densities on and in-between the stripes 
depend on Vb is shown in Fig. [3] for p = 0.774 and 
p = 0.875. Vb = corresponds to the homogeneous lat- 
tice, and the stripe and interstripe densities are equal 
there. As Vb increases, charge is driven off the stripes 
until, ultimately, for Vb > 10, the stripes are nearly 
empty. The fermions flow into the interstripe regions. 
Their density rises, going through half- filling at Vb ~ 6, 
and asymptotically increasing to a bit over unity for large 
Vb for p = 0.774. If the average density on the entire 
lattice were p = 3/4 then, for 7^ = 4, the interstripe re- 
gions would be precisely half-filled when the stripes are 
completely empty. For p = 0.875 the interstripe density 
crosses half-filling at Vb ~ 3 and approaches p ~ 1.2 as 
Vb increases. 

In Fig.|4]the spin-spin correlation c(i) is given down the 
center of the interstripe region (i.e. parallel to the stripes, 
see arrows in the inset), at p = 0.774 for Vb =4 (a) and 
Vb = 10 (b) . As seen in Fig. |3] these values correspond 
to densities slightly below and slightly above half-filling 
respectively. Despite the doping, there is a fairly robust 
antiferromagnetic pattern as T is lowered. 

In contrast, in the absence of stripes, Vb = 0, the doped 
holes are spread uniformly throughout the lattice, and for 
the same density as exhibited in Fig. |4] antiferromagnetic 
order is very short ranged, as seen in Figure [5j In the 
absence of any sort of charge inhomogeneity it would be 
very unlikely that these weak magnetic correlations could 
provide the 'pairing glue' for high temperature supercon- 
ductivity. Thus the results of Figs. |4] and [5] suggest that 
domain formation is a prerequisite for any superconduc- 
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FIG. 5: (Color online) Same as Fig. [4] except for Vb = so 
that the lattice is at homogeneous density. When there are 
no stripes the spin correlations dies out very rapidly for this 
doping of the uniform Hubbard Hamiltonian. They would be 
unlikely to be able to supply the "glue" for Cooper pairing. 

tivity which is postulated to arise from robust magnetism 
at this density. 

The presence of stripes alone is not enough to guaran- 
tee the presence of antiferromagnetic correlations. Fig. [6] 
shows the same spin-spin correlations as in Figs. |4] and 
[5] but for p = 0.875. Although spin-spin correlations are 
slightly higher at Vb = 3, where the interstripe density 
is close to one, antiferromagnetic correlations are very 
short ranged for all Vq. 

If the lattice is traversed perpendicular to the stripes, 
we expect the spin correlations to be significantly re- 
duced: the low density on the stripes does not support a 
very large moment. In Fig. m we show c(i = 2y), corre- 
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FIG. 6: (Color online) Same as Fig. [I] except for p = 0.875 
at fixed P — 5 and different values of Vq. 



sponding to a pair of sites on a line parallel to the y axis 
and traversing a stripe (see inset for a sketch). The scale 
of c(i = 2y) is roughly an order of magnitude smaller 
than c(i = 2x), as suggested should be the case by the 
preceding argument. Apart from the size of the corre- 
lation, there is another feature of crucial interest. For 
p — 0.774 and small Vq, i.e. close to the homogeneous 
limit, the spin correlation two sites away c(i — 2y) is 
positive, as expected for an antiferromagnet. However 
for Vb > 3, the sign flips and c(i = 2y) becomes nega- 
tive. This effect is strongly reduced for p = 0.875, where 
the negative values of c(i = 2y) are smaller in magnitude 
than those for p = 0.774 and only occur for Vq ^ 7. 

This 'TT-phase shift' is a promine nt exp erimental fea- 
ture of stripe physics in the cupratesSSEni. These results 
show that this shift in the sublattice order across a stripe 
is also a characteristic of the doped 2D fermion Hubbard 
model, at least in the case considered here in which the 
stripes are created through an externally imposed poten- 
tial. 

We turn now to the pairing correlations. Figure [8] 
shows the d-wave pairing correlation function along a 
stripe. As we will see below, the lowest temperatures 
achieved are well above the superconducting transition 
and therefore the superconducting correlations are short- 
ranged. Fig. [9] shows the pairing correlations between 
two neighboring sites within a stripe as a function of Vq. 
It is clear that the presence of stripes enhances pairing, 
and it is interesting to note that for both p = 0.774 and 
p — 0.875, Qpair(i) stabilizes for Vq close to the value for 
which the 'tt phase shift' takes place, namely Vq ~ 4.0 
and Vq ^ 7, respectively. 

Figure [To] shows the key result of this paper. As charge 
domains are introduced into the square lattice Hubbard 
model, the d-wave pairing vertex becomes considerably 
more attractive. Indeed, not only is its magnitude in- 
creased by a factor of three relative to the homogeneous 
system, the temperature evolution becomes increasingly 
steep. As in all existing DQMC studies of the repulsive 
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FIG. 7: (Color online) Spin correlations are shown perpen- 
dicular to the stripe, and, specifically for fixed distance i = 2y 
which crosses a stripe and varying Vb- For Vb small (the nearly 
homogeneous limit), c(i = 2y) is positive, as would arise in 
an up-down-up- ■ ■ staggered magnetic pattern. However, as 
stripes are introduced, c(i = 2y) flips sign. The magnetic 
order exhibits a 'tt phase shift' and the sublattices of the (bi- 
partite) square lattice on which the up spin electrons sit are 
reversed upon crossing a stripe. Arrows in the inset show pairs 
of sites traversing the stripe, where the spin-spin correlation 
functions are calculated. 
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FIG. 8; (Color online) D-wave pair-pair correlation functions 
along the stripes for j3 = 5, p = 0.774 and different Vq. Pairing 
correlations are shown to be short ranged. 



Hubbard Hamiltonian, the sign problem prevents access- 
ing low enough temperatures to establish a critical Tc 
where TP^ = —1, if such a superconducting transition 
does indeed occur in this model. Nevertheless, the results 
Figure [TO] are suggestive that charge domains consider- 
ably enhance the c?-wave pairing. 

We can study this same data as a function of Vq at con- 



stant temperature. Figure 11 shows results for our canon- 
ical parameters, 16x16 lattices, U — 4 and p = 0.774. We 
find that the pairing vertex becomes more and more ro- 
bust with increasing Vq. This is not completely intuitive. 
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show the d-wave pairing vertex as a function of Vo for fixed 
FIG. 9: (Color online) D-wave pair correlations on neigh- -^^^^.^^ temperature p. FPa becomes monotonically more at- 
boring sites, along the stripe, as a function of Vo- It is clear tractive. Data for 16x16 lattices (closed symbols) and 8x8 
that, as Vo is increased the d-wave pairing is enhanced along lattice (open symbols) show finite size effects are smaller than 
the stripes. statistical fiuctuations. 




FIG. 10: (Color online) The d-wave pairing vertex is shown 
as a function of temperature for U = 4, p = 0.774 and dif- 
ferent values of the externally imposed stripe potential Vo. 
For the homogeneous system, Vo = 0, and small Vo generally, 
TPd ^ 0.10. For larger Vo, TPd exceeds 0.3 in magnitude. A 
superconducting instability is signaled by TPd — > — 1. Data 
for 8x8 lattice (Vo = 4, open symbols) show finite size effects 
are negligible. 



One might expect that a maximum pairing would occur 
for Vo « 6 where, according to Fig. [3j the inter-stripe 
regions are most close to half-filling and hence antiferro- 
magnetic correlations are most strong. Alternately, as Vq 
becomes very large the particle density within the stripes 
drops to zero. If the physical picture is that of pairing of 
mobile carriers in the doped region driven by spin corre- 
lations in the half-filled domains, at large Vq the density 
of these carriers becomes small, and one would again ex- 
pect TPd to turn over. 

Recent DCA calculations^ indeed reveal an initial en- 
hancement of pairing with the introduction of stripes. 



followed by a fall-off as the above physical arguments sug- 
gest. However, this non-monotonicity is observed only at 
very low temperatures T ^ 0.05 quite close to Tc. For 
higher T, closer to the range studied here, there is no 
sign of the d-wave eigenvalue coming down at large Vq. 
Although the lower temperature scales is the most likely 
explanation for the DCA non-monotonicity, it is also pos- 
sible that, in the work presented here, the decline of pair- 
ing correlations in the doped stripes is compensated by 
an increase in the interstripe domains, which are shifted 
away from half- filling at large Vq. It seems clear that for 
large Vq and density p = 0.750 where the interstripe re- 
gions are precisely half filled, the pairing signal would be 
forced to be small, and hence that a turnover such as is 
seen in DCA calculations should occur. 

The physics of the half-filled homogeneous Hubbard 
model on a square lattice is believed not to be highly 
sensitive to the interaction strength U. That is, the 
ground state is an antiferromagnetic insulator for all U, 
although the precise nature of the phase evolves from 
a weak coupling regime described by spin density wave 
physics to a strong coupling Mott insulator. In order to 
assess whether the enhancement of d wave pairing due 
to striped formation is similarly generic to different U or 
specific to U = 4, we show data for TPd as a function 
of temperature T in Fig. [12] for U ^ 6. The same basic 
evolution is observed as in Fig. |10[ with the vertex being 
enhanced by Vq both in magnitude and also in the steep- 
ness of its evolution as T is lowered. A comparison of the 
data ranges of Figs. [lO] and [12] also reveals some of the 
limitations of DQMC. As U gets larger the sign problem 
grows, as do the fluctuations (error bars). For U — 6 the 
lowest accessible temperature is T ^ 0.3 {(3 — 3), com- 
pared to T ~ 0.2 (/3 = 5) at [/ = 4. Smaller interaction 
strengths, e.g. C/ ~ 2, can readily be simulated, but tend 
to harbor large finite size effects which are the remnants 
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FIG. 12: (Color online) The pairing vertex is shown as a 
function of temperature for the same parameters as Fig. |10| 
except for a larger on-site repulsion U = 6 rather than U = 
4, and 8x8 lattices. The generic behavior is similar: charge 
inhomogeneities make the d-wave vertex more attractive. 




FIG. 13: (Color online) Pairing vertex for the V = 2 case, 
where the stripes are separated by a single interstripe row. 
Unlike the "P = 4 case where a three-site-wide interstripe 
region separates the stripes, increasing charge inhomogeneity 
is detrimental to pairing. In fact, the d-wave vertex becomes 
repulsive for Vb > 2. 



of high degeneracies in the noninteracting energy levels 
on a square lattice. 

We next explore a different lattice periodicity, V = 2. 
Unlike the case of different interaction strengths, where 
the qualitative physics remains unchanged, V = 2 be- 
haves in a very different manner from 7^ = 4, as seen in 
Fig. |13| In this case, increasing Vq reduces the magni- 
tude of TPd, and by the time Vg = 3 the vertex has even 
changed sign and become repulsive. It should be noted 
that because for V — 2 there are fewer stripes which do 
not feel the Vq potential to absorb the fermions as Vo in- 
creases, their density is increased above p = 1, where AF 
correlations are most evident, far more easily than in the 

= 4 case shown in Fig. |3] This is the probable cause 
for the decrease in d wave pairing correlations. 



FIG. 14: (Color online) For a total density corresponding 
to close to the 'optimal' doping of the cuprate superconduc- 
tors, p = 0.875, the pairing vertex, while remaining attractive, 
shows somewhat less enhancement as Vb is turned on. 



Our final check of the robustness of the enhancement 
of pairing by stripes is to explore a different density 
p = 0.875, which is optimal for pairing in the absence of 



stripes. Fig. 14 indicates that the temperature evolution 
for non-zero Vq is essentially unchanged from the homo- 
geneous case. Thus, even though at Vb = the product 
TPd is larger for p — 0.875 than for p = 0.774, this is no 
longer the case for nonzero Vq. Indeed, stripes can ulti- 
mately make the pairing vertex substantially larger for 
the lower density. The strong antiferromagnetic correla- 
tions induced by the stripes at V = A provide the "glue" 
for pairing a.t p = 0.774, whereas at p = 0.875 this "glue" 
is not present. 

We have focussed thus far on the density, spin, and 
pairing correlations. F ig . 1 1 5 1 examines the total energy as 
a function of Vq on 16x16 lattices. Although we have im- 
posed Vq in our Hamiltonian, this computation of the en- 
ergy provides a crude measure of the tendency for spon- 
taneous stripe formation. The monotonically decreasing 
behavior of the energy with Vq, suggests that a maxi- 
mization of charge imbalance is favored. When we study 
the effect of the periodicity of the stripes, as shown in 
Fig. [16] on a 12x12 lattice, we see that 7^ = 3,4 have 
much lower energies, suggesting their formation might 
be favored. These periodicities have densities on the lines 
which are not subject to the additional potential Vq rela- 
tively close to half-filling, and hence they have the largest 
antiferromagnetic correlations. 



IV. CONCLUSIONS 

The external imposition of stripes via the introduc- 
tion of a linear pattern of reduced chemical potential has 
been shown to result in a significant enhancement of the 
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FIG. 15: (Color online) The total energy as a function of Vo 
for 16x16 lattices with (7 = 4, /3 = 5, P = 4, and p = 0.875 
(squares) and p = 0.774 (circles). 




FIG. 16: (Color online) The total energy for different dis- 
tances between stripes. The minimum for "P = 3, 4 is associ- 
ated with the fact that the density of the Vb = rows is close 
to half-filling for total density p — 0.774. 



d-wave pairing vertex of the two dimensional Hubbard 
Hamiltonian. When the overall density and periodicity 
of the stripes are such that the density in the interstripe 
region is close to one, antiferromagnetic correlations are 
also made larger, and exhibit a 7r-phase shift across the 
stripes. Both the 7r-phase shift and the growth in the su- 
perconducting response occur only when the charge inho- 
mogeneity is sufficiently large, specifically when the ad- 



ditional inhomogeneous site potential Vq exceeds roughly 
three times the hopping. 

This enhancement of superconductivity has previously 
been observed in the closely related dynamical cluster ap- 
proximation treatment of the two-dimensional Hubbard 
Hamiltonian, again in the case when a site potential Vq 
was introduced externally. In this situation, the modu- 
lation was chosen to be broader than the purely one di- 
mensional pattern explored here. The observation of an 
optimal stripe potential in the DCA calculations might 
be associated with this difference. It is also interesting to 
note the possible differences betwe en pin ned and fluctu- 
ating stripes. There are suggestion^--'^^"^^' that the mo- 
tion of charge/spin domain walls is important to the en- 
hancement of superconductivity, whereas frozen stripes, 
such as created by Nd doping, is inimical to pairing. The 
studies in the present paper, as well as the earlier DCA 
work of Maier eta¥^ indicate that the situation may not 
be quite so straightforward and that, in fact, the stripes 
produced by an external potential might also be able to 
enhance superconductivity. 

We have motivated our form for the externally imposed 
stripe potential at T' = 4 as producing the charge/spin 
patterns suggested by neutron scatterin^^ll and DMRG 
calculations—' Our pattern contains equal Fourier com- 
ponents for all wave vectors Q„ = 2TTn/V. In the 
DCA worlPSl the effect of different Fourier components 
Q = tt/2 and Q = tt/A was explicitly isolated, with the 
former showing little effect on pairing and the latter driv- 
ing significant enhancement. 

Spontaneous stripe formation in the doped Hubbard 
model, if it occurs, takes place at temperatures below 
those accessible to DQMC simulations, which are lim- 
ited by the sign problem'^ to temperatures greater than 
roughly one-fortieth of the noninteracting bandwidth. 
It would be interesting also to explore the possible en- 
hancement of pairing at these temperature scales by 
other types of charge inhomogencities such as checker- 
board patter nj^ and in the presence of nonmagnetic 
disordeiE^. 
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